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Abstract 

We present a calculation of the branching ratio of orthopositronium decay into an 
invisible mode, which is done in the context of Mirror World models, where ordinary 
positronium can disappear from our world due to oscillation into its mirror twin. In 
this revision we clarify some formulas and approximations used previously, correct 
them at some places, add new effects relevant for a feasible experiment and finally 
perform a combined analysis. We include into consideration various effects due to 
external magnetic and electric fields, collisions with cavity walls and scattering off gas 
atoms in the cavity. Oscillations of the Rydberg positroniums are also considered. To 
perform a numerical estimates in a realistic case we wrote computer code, which can 
be adopted in any experimental setup. Its work is illustrated with an example of a 
planned positronium experiment within the AEgIS project. 

1 Introduction 

Since discovery of parity violation in particle physics in 1957 [1] its origin in Nature remains 
a mystery. An idea for how to retain the mirror symmetry based on new particles of opposite 
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chiralities was put forward long ago [2] and developed to the idea of two coexisting parallel 
Worlds in [3] (see also [4]), where mirror transformation can be realized as the composition 
of CP and particle-H-mirror particle transformations. 

The idea of a Mirror World has a long and interesting history nicely recounted in the 
review [5]. Today the Mirror World model is one of the most attractive possible extensions 
of the Standard Model of particle physics. This theory is based on the assumption that 
in addition to the Standard Model its mirror image exists and corresponding new particles 
interact very weakly with SM fields. Thus, the theory is based on the gauge group G x G", 
where G = SU{3) x SU{2)l x f/(l)y, G' = SU{3y x 5f/(2)'R x f/(l);.. As it was pointed 
out before, originally the deep theoretical motivation for this extension was originally a 
restoration of the mirror symmetry: P-transformation changes ordinary left particles to 
mirror right ones [2]. If this symmetry is exact then every particle has its mirror partner 
(twin) with the same mass and lifetime but is charged under different gauge groups. 

Certainly, if mirror particles exist they inevitably interact with ordinary particles gravita- 
tionally. But in general these models allow for two kinds of direct renormalizable interactions 
— portals — between the Standard Model (SM) and the Mirror World particles. One of 
them, called the Higgs portal [6], can be probed in high energy collider experiments (see e.g. 
Ref. [7]). Here we discuss implications of the second type of interactions between SM and 
Mirror World particles, dubbed as the Abelian portal, 

2 COS 9w COS e',^ ^" ' ^ ^ 

where S^,^ and 9w are the U{1)y field strength tensor and the weak mixing angle, respec- 
tively, and B'^j^ and 9{y are the same quantities in the mirror sector. Studying of this 
interaction calls for another sort of experiment. 

The most sensitive way to search for the Mirror World in a laboratory is looking for the 
disappearance of ordinary orthopositronium (oPs) resulting from oscillation to the mirror 
twin and its subsequent decay to mirror photons. The oscillation becomes possible due to 
a photon - mirror photon mixing hidden in Eq. (1). Similar processes for bound states of 
quarks (vector mesons) are much less sensitive to the Abelian portal (1). 

A new experiment [8] is proposed to search for the oPs disappearance with the sensitivity 
to branching ratio at the level of 10~^ — 10~^. This refers to the value of the mixing parameter 
entering (1) e ~ lO"^. Quite remarkebly, it is the first time when experimentalists have a 
chance to enter a region of parameter space which is both phenomenologically viable and 
interesting for cosmological applications. 
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Indeed, the mirror matter can play very important role in cosmology and astrophysics 
(see, e.g. Refs [9, 10, 11]). Mirror matter can behave like dark matter on the astrophysical 
and cosmological scales. With Mirror symmetry (slightly) violated (say, by different vacuum 
expectation values of our Higgs field and Mirror Higgs field) the phenomenology becomes 
more fascinating. Indeed, photon-paraphoton mixing from the Abelian portal (1) implies 
the presence of millicharged particles: mirror particles effectively carry a tiny charge with 
respect to electromagnetism. Thus, mirror matter particles can be directly produced in col- 
lider experiments [12], be tested with electromagnetic precision measurements [13], change 
the CMB anisotropy picture due to participation in primordial plasma dynamics at recom- 
bination [14], contribute to supernovae explosions [15]. The idea of Mirror World naturally 
explains the presence of light sterile neutrino(s) at the eV-scale, as suggested by neutrino 
anomaly announced by the LSND experiment [16] and is now in agreement with the com- 
bined analysis of cosmological data [17, 18, 19]. Indeed, since we have sub-eV neutrinos in 
our sector, which within the SM acquire masses from dimension-5 interaction terms, then 
mirror (sterile) neutrino masses are of the same order. Mirror symmetry provides similar 
dimension-5 operators which give masses to mirror neutrinos and may yield mixing with 
neutrinos from our sector. 

In this paper we revise positronium oscillations into its mirror twin. The paper is orga- 
nized as follows. We start in Section 2 with the discussion of positronium oscillations in a 
vacuum. Further, in Section 3 we describe the relevant parameters of the proposed experi- 
ment [8] (Section 3.1), present the combined analysis of positronium scattering off the gas 
atoms (Section 3.2) and the cavity walls (Section 3.3), and take into account the interaction 
with electric and magnetic fields (Section 3.4). Numerical estimates are given for a particular 
setup of the proposed experiment [8]. Section 4 is devoted to the study of oscillations of the 
Rydberg (highly excited) positroniums. Section 5 contains conclusions and a discussion of 
open problems. 

2 Warming-up task: orthopositronium oscillations in 
vacuum 

We begin with a simplified picture of vacuum oscillations. Hereafter if not stated otherwise 
we consider a model with the Mirror World, where the mirror symmetry is exact: values 
of all mirror masses and mirror coupling constants coincide with those of the corresponding 
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originals. Positronium couples to its mirror twin due to photon-paraphoton mixing emerging 
from the Abelian t/(l)y-portal (1), 

Unt = \ F'^^F^^ , (2) 

where F^,^ and F'^^ are electromagnetic and paraelectromagnetic field strength tensors, re- 
spectively. Orthopositronium (oPs) can oscillate to its twin oPs' via virtual photon mixed 
with virtual paraphoton due to coupling (2). The corresponding transition matrix element 
reads 

r 

(oPs I Cnt I oPs') = 27re/ = - , (3) 

where parameter 

f^S.7x 10^ MHz 

is determined [20] by the one-photon annihilation diagram involving orthopositronium. Ex- 
perimental constraints on mixing parameter e have been recently reviewed in Ref. [10]. The 
first experimental limits on e from oPs physics were discussed in Refs. [21, 22]. Searches 
for an invisible decay of orthopositronium have a rather long history [23, 24, 25]. Nonob- 
servation of positronium disappearance oPs invisible places a direct limit on the mixing 
parameter [26] 

e < 1.55 X 10"^ . 

The strongest indirect limit on e in a model with exact mirror symmetry comes from 
cosmology [10], 

e < 3 X 10-^° . (4) 

It ensures that our and Mirror Worlds have never come to equilibrium in the early Universe; 
then with somewhat colder mirror plasma the mirror baryons serve as cold dark matter in 
the Universe. Remarkably, the limit (4) is not far from the estimate [27] 

~ (7 ± 3) X 10-^° , (5) 

required to explain the annual modulation signal, in direct dark matter searches performed 
by CoGeNT and DAMA experiments. Here C,o' is a mass fraction of the mirror oxygen in 
dark matter particles. Strictly speaking, there is no overlap between the two ranges (5) and 
(4). The gap grows on account of the fact that observation of the matter distribution in 
the Bullet cluster (IE 0657-558) can be reconciled with the Mirror dark matter only if the 
largest fraction of dark matter particles are confined in compact macroscopic objects (mirror 
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asteroids, stars, etc) [28]. Then, the dark matter particle flux should be lower than what 
was supposed in [27] and hence the required DAMA/CoGeNT signal interval of the mixing 
parameter is higher than (5). At the same time, the analyses of new data on primordial 
Nucleosynthesis [29] and cosmic microwave background anisotropy [30] suggest a somewhat 
higher rate of the Universe expansion which can be explained by a contribution of mirror 
matter. This relaxes the cosmological limit (4) to some extent. 

In the nonsymmetric model, the cosmological bound is generally much weaker. For 
instance, let the only source of asymmetry be different values of the Higgs boson vacuum 
expectation values, (H) ^ {H'). Then in the hierarchical case (H) <C {H') one obtains [10] 

,<3xlO-'x^, 

from the succesfulness of the Big Bang Nucleosynthesis. 

To summarize the limits discussed above, in what follows we adopt 

e = 1 X lO-'^ (6) 

as a reference number for numerical estimates. This particular value is of considerable 
interest, since it has been argued [31] that the DAMA/LIBRA periodic signal [32] and obser- 
vations by CoGeNT [33] can be explained within the Mirror World model with a mixing of 
the order (6) by galactic mirror baryons, while other relevant experiments [34] were blind to 
this dark matter. In any case, if mirror matter plays the role of dark matter particles, the di- 
rect limit on mirror mixing should be of order (6). Moreover, it has been argued that studies 
of positronium in the proposed experiment [8] will be sensitive right up to the values of mix- 
ing parameter. Thus, for the flrst time the experiment enters into the phenomenologically 
interesting and cosmologically allowed region of the model parameter space. 

For the reference mixing (6) the value of matrix element responsible for oPs ■H- oPs' 
oscillations (3) is 

6 « 1090 s"^ ^ 7.3 X 10"^^ eV , (7) 

which is much lower than the rate of orthopositronium decay (mostly decay to three photons) 
[35, 36] (see Ref .[37] for more precise and recent theoretical calculations), 

7 ^ 7.040 X 10^ s"^ ^ 4.634 x 10"^ eV . (8) 

Hereafter we work in the "flavor" basis: the evolution of the system oPs-oPs' in a vac- 
uum can be described by the Schrodinger equation for the doublet wave function — 
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T 

(V'oPs,V'oPs') , 

= (9) 

where the Hamiltonian in the positronium rest frame accounts for the positronium and mirror 
positronium decays and their mutual oscillations, 

I, (10) 

where E is the kinetic energy of oPs and oPs'. When solving Equation (9) for the wave 
function one obtains the probability to have mirror positronium instead of initial positronium 
by the elapsed time t, 

P (oPs ^ oPs') = e-^* sin^ y . (11) 

Mirror orthopositronium decays mainly into three mirror photons and the probability of oPs 
disappearance is given by [38] 

/•oo p ^ p 

Br (oPs invisible) = 7 / P (oPs oPs') dt = ^ — , (12) 

7o 27^ + 2Y 

since 5 < 7, cf. Eqs. (7) and (8). 



3 Realistic consideration 
3.1 Realistic experimental setup 

We proceed with the study of complications in the positronium description arisen in a real 
experiment. Indeed, the realistic case of positronium oscillations is much more involved. In 
our analysis we consider the general setup of experiments proposed in Ref. [39] and then 
further developed in Ref. [8] with more reliable estimates of the sensitivity to the mixing 
parameter. 

In the experiment [8] oPs states will be formed in a thin nanoporous Si02 target placed 
on the bottom of the vacuum cavity (10 cm in diameter and 10 cm in height). After the 
formation a significant fraction of the oPs can escape inside the cavity and can become 
almost completely thermalized with the average temperature equal to the temperature of the 
target [35] . The cavity is surrounded by an almost hermetic 47r electromagnetic calorimeter 
to detect annihilation photons. The experimental signature of the oPs-oPs' oscillations is 
the absence of energy deposition in the calorimeter expected from the ordinary positronium 



6 



decays. The typical residual gas pressure is above lO^^Torr. There is an optional external 
electric field E < lOOkV/cm and magnetic field B < 10'^ G (exceptionally up to < 10^ G). 

Note that the experiment [8] can also be preformed, e.g., in the framework of the AEgIS 
project (Antihydrogen Experiment: Gravity, Interferometry, Spectroscopy) at CERN [40, 
41]. The main goal of AEgIS experiment is a mass production of antihydrogen. In this ex- 
periment orthopositroniums are a source of slow positrons that recombined with antiprotons. 
There is a proposed technology of adopting Rydberg (highly excited) positroniums instead 
of the usual ground-level orthopositroniums: the former live longer and positroniums move 
slower inside. Thus, in the AEgIS experiment, not only ground states will by confined within 
the cavity, but also excited oPs with the Rydberg numbers n < 30. 



3.2 Oscillations in gas 

Let us consider the role of positronium interaction with gas. Generally, the time-evolution 
of the system oPs-oPs' is described by density matrix 

p{t)^ fd'x^^^^ f£x( ] . (13) 

J \ V'oPsV'oPs' V'oPs'V'oPs' / 

The density matrix (13) solves the following equation [38] 

^ = -mp + ipn^ + 2nnvjd cos 9 F{e) p F\9) , (14) 
with the Hamiltonian, cf. (10), 

5/2 E'-^i/2)' ^''^ 

where E and E' are the positronium and mirror positronium energies, 7 and 7' are their 
widths^, k refers to the value of the positronium 3- momentum, and v is the mean relative 
velocity between the positronium and gas molecules. The positronium scattering off gas is 
described by the matrix 



^Different notations for oPs and oPs' are used {E, E' for energies and 7,7' for widths) because corre- 
sponding quantities can differ from one another in external electromagnetic field and for the case of broken 
mirror symmetry. 
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where is a scattering angle; in the cavity mirror positronium does not scatter off ordinary 
matter given the smaU mixing (6). The first term in the Hamiltonian (15) accounts for the 
presence of gas in the system, which changes the energy of the positronium state. There 
n stands for the number density of gas particles, so that the frequency of the positronium 
scatterings off gas w is determined by the cross section a or (with help of the optical theorem) 
by the forward scattering amplitude /(O) as 

47r 

w = anv — — nv lmf{0) . (17) 

In what follows, we resort to the case of slow scatterings, where the scattering amplitude 
does not actually depend on the scattering angle 9: i.e. scattering is saturated in s-wave. 
The real part of the scattering amplitude is positive^ and by analogy with Eq. (17), can be 
conveniently parametrized as 

An 

«;fle = ynT;Re/(0) . (18) 

In this case the problem is solved analytically as follows. Once positronium is produced (at 
t — 0), the system is in the pure flavor state, 

^(0)=^^j , and hence p(0)=| ] . (19) 

We parametrize the density matrix at arbitrary moment as 

P-l I . (20) 




x-iy p2 



then Eq. (14) yields 



^ = -7Pi - 5y (21a) 

^ = 5y-7'p2 (21b) 

| = ^(pi-P2)+(^-A)a;-i(7 + y + ^)y (21c) 

^ = -1(7 + 7' + ^)^+ (A-^)y (21d) 



where A = E' - E. 



^At large distances r positronium is attracted by an atom, i.e., it propagates in the negative potential 
U (r) < 0, hence in the Born approximation for nonrelativistic particles we obtain positive value for the 



amplitude (see, e.g. [42]), / w -2m U{r)r'^dr > 0. 
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We are interested in the total probability for the orthopositronium to disappear, which 
under study reads 

Br (oPs ^ invisible) = ,00 {° ^'fl ■ (22) 

Jo Pidt + Jo P2dt 

One can integrate the system (21), with the initial condition (19), and introducing notations 

/•oo roo 1*00 

Y= ydt , X= xdt , Pi= pidt , i = 1, 2 , 
Jo Jo Jo 

arrive at a linear system of algebraic equations: 

-7P1 -6Y = -1 
SY - = 

Kt^ n \ , f^Re , \ ^ If , , , V „ (23) 
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(Pi-P2)+(^-A)x-1(7 + 7' + ^)>^ = 
1(7 + 7' + ^)^+(A-^)f = 



2 

Then the oscillation probability (or probability to disappear) (22) is given by 

Br (oPs invisible) = . (24) 

Pi + P2 

Let us introduce the following notations: 

r = ^(7 + 7' + w), An, = A-WRe/'2. (25) 

Taking into account the hierarchy 5 <^ F, the solution of the system (23) can be simplified 
and we obtain P^ + P2 = ^ (l + O (^^^^ and for the probabihty [38], 

Br (oPs ^ invisible) = y ^ ^ ^ (f^) ) " ^^^^ 

Note that in the "flavor" basis, the case of an exact mirror symmetry implies equal energies 
and widths of positronium and its twin in the absence of external electromagnetic fields, 
hence A = and 7' = 7. 

For the numerical estimate, we apply the obtained formula to a particular case of the 
experiment described in Section 3.1. Here the cavity is proposed to be filled with nitrogen 
N2 at room temperature Tq, so that 

A;bTo = 0.025 eV . (27) 
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The momentum transfer cross section of positronium on nitrogen was measured in Ref . [43] 
as follows: 

/I 
(1 - cos 6) -^dQ = (35 ± 8) X lO'^^cm^ , (28) 
ail 

and remains almost constant with temperature varying at ksTo < 0.3 eV. And if scattering 
amplitude does not depend on scattering angle we obtain for the total cross section after 
integrating over the angle variables a = J ^dQ In accordance with the feasible 

experimental setup [8] we assume that orthopositroniums are emitted from the target with 
a temperature of T^pg ~ Tq, close to the nitrogen temperature (27). Thus, the positronium 
average momentum is A; = Jtt.oPsI'oPs where positronium average velocity equals 

VoP.= \ = 8.1x10*' cms X W . 29 

V moPs V 0.025 eV 

In principle, there may be processes of pick-off annihilation and ortho-to-parapositronium 
conversion. But the corresponding probabilities are neghgibly small (see, e.g. Ref. [44]) 
at the positronium energies and the gas pressures relevant for the class of experiments in 
question [8]. 

In what follows, we set ToPs = Tq. One observes, that VoPs is much higher than the 
velocity of nitrogen atoms, so one can replace v in Eq. (17) with v^ps in order to estimate 
the collision frequency and, in particular, to check that for the cross section (28), the slow 
scattering approximation is valid indeed and the scattering amplitude does not depend on 
the scattering angle as we have assumed. For corresponding energies of colliding particles 
the cross section is mostly saturated by the real part of the amplitude, so 

|Re/(0)| ^ y/a/An . 

The nitrogen pressure Pjv2 is related to the nitrogen number density riiy.^ by P^^ = 
n]\f^kBTN2- Obviously, the presence of gas suppresses oscillations (see Eq. (26)), if its density 
is high enough for positronium to scatter off atoms at least once. This effect can be under- 
stood as the loss of coherence due to a different angular distribution of the positronium and 
its mirror twin after collision, which is encoded in Eq. (16) (see also Ref. [45] for discussion 
of similar effects). For the reference temperature (27) this implies a pressure Pn^ of order 
10~^ Torr or higher. We put all the formulas above into Eq. (26) with the reference value of 
the mixing parameter (6). When introducing dimensionless variables V = PN2/ (10~^ Torr) 
and T = ks Tj^^/ (0.025 eV) we write for the invisible branching ratio 

Br (oPs invisible) = 1.2 x 10- x (^)\ ^ ^Ml^vn^A^^^^^^ x V /T)^ " 

(30) 



10 



This dependence can be used [21] to check against the possible systematics once the signal 
is found. 



3.3 Oscillations in a finite volume 

In a laboratory orthopositronium is trapped in a cavity. Then if positronium velocities are 
high enough to reach the walls before decay, the interaction with wall material modifies the 
(possible) oscillations into mirror positronium. The reason is simple: positronium interacts 
with walls, while its twin does not and if produced in a cavity flies away freely. 

Let us use the instant approximation for the positronium interaction with a wall. Then, 
between the scatterings off the wall, the oscillations proceed as in an infinite volume case, but 
when positronium hits the wall, nondiagonal elements of the density matrix (13) nullify. The 
oscillations stop at this point, since the mirror positronium flies away, while the positronium 
gets reflected, and hence their wave functions get separated in space, and the coherence gets 
lost. This is the main observation. 

To account for this instant decoherence properly, let us consider the oscillations of a 
stable positronium. Then, the evolution of the density matrix elements is described by the 
following equations: 

IF = P^'') 

(310) 

and X remains zero with initial condition (19), which is the pure positronium state. 

It is convenient to introduce the variable s = pi — p2 and parametrize the time interval 
between i — 1-th and i-th reflections as Tf, the system starts to evolve at t = being in the 
positronium state (19). Then before the flrst reflection the solution of the system (31) is 

t < Ti : s = cos to , y = - smto . 

At i = Ti we nullify y (ti) = 0, but the diagonal variable remains intact, s{t\) — costi(5. 
Then, the evolution of the system between the flrst and second reflections proceeds as follows, 

Ti < i < Ti + T2 : s — cos Ti5 COS (i — Ti)(5 , y — ^'^^ '^'^^ — ■ 
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It is straightforward to see that right after the n-th reflection the solution reads 

n 

5 (T"n) = Yi ^1^^ ' yi^n) = ■ 



k=l 

When introducing the average time between reflections 

t 

n 

one finds that the positronium velocity (29) and the typical size L ~ 10 cm of the detector 
volume are related as (r) ~ L/voPs- For the reference value of the oscillation rate (7), the 
average reflection rate obeys 

(r)5 < 1 . 

Then we can use the exponential approximation to obtain the smooth solution 

Thus, collisions with walls result in the exponential suppression of oscillations. 

One can mimic the suppression (32) by introducing an additional suppression to the 
Equation (31c), 

dy S 

Then, at 5/w <C 1 one has 

/ 5H 

s = exp 

\ w 

and matching (32) is achieved with 



w=|^. (34) 

Thus, comparing (33) to (21c) one accounts for the collisions with walls by replacing (see 
(25)) 

r ^ r + w 

in the formula for oscillations (26). 

To support our approach, let us consider a limiting case when the time intervals between 
the coUisions with the walls are all equal, t^+i — = t, then (r^)^/^ — {t) — t and one 
obtains for the oscillation probability 
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1 / (5^ Tt \ 1 . . 

2 (l-e-^j (35) 
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This result is in agreement with Rcfs. [46, 47] for the neutron-antineutron and neutron-mirror 
neutron systems where any neutron collides for many times. On the contrary, the obtained 
result deviates from a similar one in Ref. [8] where w = l/(r) instead of (34), hence it is 
two times smaller. The origin of this term in the corresponding equations of Ref. [8] was not 
explained in that paper. One can obtain this answer with our method provided the random 
variable r is distributed in accordance with the Poisson statistics, which we find unjustified 
when the positronium trajectories are reasonably tractable. 

However, in general case, the wall collision rate is not a quantity which can be well- 
defined in any real experiment. For example, the authors of Ref. [8] argue that the number 
of collisions in the proposed experiment is of order unity, i.e. ~ 1. In this case, the 
gas pressure, the geometry of cavity, and the properties of the oPs source (like energy and 
angle distributions) become important and the wall collision rate becomes a largely uncertain 
quantity (e.g., it depends on the direction of positronium propagation, etc.). 

To take into account all those effects and make accurate predictions for the probability 
of positronium disappearance, a numerical calculation should be performed. We apply the 
Monte-Carlo approach to simulate the evolution of the oPs-oPs' system in the following way. 
The density matrix is calculated numerically by solving Eq. (21) with initial conditions (19). 
We use the geometry of the positronium cavity described in Section 3.1. In our numerical 
simulations, we suppose that the positroniums in the initial state have isotropic nonrela- 
tivistic Maxwell velocity distribution with a temperature of TqPs (equal to the temperature 
of gas), which is reasonably close to the realistic distribution observed in Ref. [35]. When 
collision with the wall takes place we change the density matrix in the following way 

\p2i p22 J \ y 

which corresponds to total disappearance of mirror positronium from the cavity. The deco- 
herence due to gas coUisions is already taken into account in Eq. (21). However, coUisions 
with gas also change the direction of positronium motion. So, we simulate time intervals 
between the collisions according to its average value (17): the probability of the positronium 
experiencing an interaction in small time interval dt equals Pint — Pu (l — e"""'*). After the 
collision with a gas atom, the direction of the positronium velocity changes according to the 
isotropic distribution. We neglect small changes in the absolute velocity in the scattering off 
gas atoms. To get feehng of the influence of wall coUision let us consider the same topology 
of experiment cavity as it is described in Section 3.1 but with different absolute size, i.e. 
we take vacuum cavity of L cm in diameter and L cm in height. In Fig. 1 we present the 
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Figure 1: The dependence of the branching ratio on the characteristic size of experiment 
cavity L at different values of gas pressure. We chose e = 10~^, ksToPs = 0.025 eV. 

branching ratio dependence on the characteristic size L of cavity where orthopositroniums 
decay. 

The dependencies of the branching ratio of orthopositronium disappearance on gas tem- 
perature and gas pressure are presented in Fig. 2. To see clearly the effect of wall collisions 
we also plot the pressure dependencies obtained with Eq. (30) for comparison. Obviously, 
the cavity finite size effect as well as the gas effect are always in a reduction of positronium 
disappearance branching ratio: collisions prevent positronium from oscillations to its twin. 

3.4 Oscillations in static electric and magnetic fields: Zeeman and 
Stark effects 

It is well known that in an external static magnetic field B the state of orthopositronium 
with 171 = mixes with parapositronium [48, 49] and these states show Zeeman shifts in 
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Figure 2: The dependence of the branching ratio of the positronium disappearance on gas 
(here is nitrogen) pressure (left panel) and gas temperature (right panel) with (points) and 
without (solid lines, Eq. (30)) wall collisions. The positronium temperature equals nitrogen 
temperature and we set e = 10~^. Note that Br(oPs — invisible) scales as e^. 

their energies and decay rates^ 



A 



± 



-A 



HFS 



± 



A2 



+ (4/ioS) 



(37) 



± 



^ I 7 + 7p ± 



(7 - 7p) ^hfs 



^"hfs 



(38) 



.1 X 10^ s"^ is 



(4/ioi?)2 

where Ahfs ~ 8-4 x 10~^ eV is hyperfine energy splitting [50], 7p 
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its contribution to oPs-oPs' oscillations is negligibly small. When the magnetic field is weak, 
4/xo-B ^ Ahfs, the shifts in the energy levels and decay rates read 



A, 



A 



HFS 



Ahfs 



7 + 



Ahfs 



A_ 



7p , 



-A 



HFS 



2fioB 
Ahfs 



Ahfs 
2 



7 + 



2fioB 
Ahfs 



If ■ 



•^In the limit i? — > the state "+" goes to the orthopositroniuni state while "— " to parapositronium. 
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Thus, when the weak magnetic field is apphed the oPs state becomes a superposition of 
"+" and "— " states of which "— " decays quickly, while "+" is approximately m = state 
of the orthopositronium which acquires shifts in its energy and decay rate. Numerically we 
estimate ^ 

A+ ^ 2.5 X 10^ s"^ X I -r^ ) . (39) 



.100 

The mixing between new state and the mirror twin of m = state of the orthopositronium 
also gets corrections of order [hqB / ^hfsY which can be neglected. 

There is also the quadratic Zeeman effect in positronium [51]. The corresponding energy 
shift 

l^u^2a'alB'^lQU-^x(^-^^^ , (40) 

where ao — 0.1 nm is the Bohr radius, applies equally to all four spin states of the lowest 
energy and can be considerable for large magnetic field. 

In an external static electric field E all three orthopositronium states get shifted equally 
due to the Stark effect [42], 

A. = -\a„e„E^ = -1.85 X lO-eV x (^5^3^)^ "2.8 • 10= s"' x (^5^3^)' , 

(41) 

where ckq = |47raQ determines polarizability of orthopositronium, and eo refers to the free 
space permittivity. With a reasonable field strength available in a laboratory (in particular, 
in the AEgIS experiment), E ~ lOOkV/cm, the energy splitting due to the Stark effect (41) 
can be considerable. 

Let us summarize the effect of the external electromagnetic field. For large magnetic field, 
m = state immediately decays and does not contribute to the oscillations. Hence, in case 
of unpolarized orthopositroniums this results in decreasing of total oscillation probability by 
a factor 2/3. This behavior gives a tool to check for possible systematics in the experiment, 
where magnetic fields as large as 10 T are available (as, in particular, in the experiment [8]), 
which is quite useful especially if evidence for the oscillations is found. In the weak magnetic 
and electric fields the m = oPs state acquires total energy shift 

Ao = A+ + Ac + A5 , (42) 

and its decay width becomes 70 = r+. The oPs states with m = ±1 get energy shift 

Ai = Ao + A5, (43) 
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while their decay rates remain intact. 

Averaging the disappearance probabihty over states with different angular momentum 
components 

Br oPs ^ invisible = , , xr— , 44 

one obtains 

Br (oPs ^ invisible) = ^ ( ^ , ° „ + 2^-^ — — , (45) 

where Pq = |(P+ + 7 + w). 

The results for invisible branching ratio of the orthopositronium in the presence of electric 
and magnetic fields with different values of gas pressure are presented in Figures 3 and 4, 
respectively. Here we also include the effect of the wall collisions in the geometry described 




le-12 I ' ' ' ' 1 

20 40 60 80 100 

E, kV/cm 

Figure 3: The dependence of the branching ratio of the positronium disappearance on the 
external electric field E at different values of gas pressure. We chose e = 10~^, fc^ToPs = 
0.025 eV and the magnetic field 5 = 0. 

in Section 3.1. We see that generically presence of an external electromagnetic field leads to 
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le-08 




B, Gauss 

Figure 4: The dependence of the branching ratio of the positronium disappearance on the 
magnetic field B at different values of gas pressure. We chose e = 10~^, fc^ToPs = 0.025 eV 
and the electric field E = 0. 

decreasing of the oscillation rate and thus to deterioration of the sensitivity of experiment 
to the Mirror World physics. Some nontrivial dependence on the magnetic field stems from 
the fact that both Zeeman effect contributions to the energy splitting have the sign which is 
opposite to the contribution of the coherent scattering (15). As we have already mentioned 
the dependence on the external parameters like electric or magnetic fields can be used to get 
rid of systematic uncertainties or (in case of positive signal) to get evidence that the observed 
effect is due to the oscillation nature of positronium disappearance rather than direct decay 
into some invisible particles. 

Let us note in passing, that in a model with broken mirror symmetry where the positron- 
ium and its mirror twin obtain mass splitting larger than orthopositronium width, oscillations 
are absent because of the loss of coherence, and the branching ratio of an invisible mode is 
strongly suppressed (see Eq. (26)), 

S'^T 1 

Br (oPs — )■ invisible) ~ — — . 

2 7 

However, with magnetic field tuned to cancel this mass splitting by using the Zeeman effect. 
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one still has a possibility of the Mirror World hunting via the positronium portal. 



4 Rydberg states and models with vacuum mass split- 
ting between twins 

Let us proceed with discussion of oscillations of (highly) excited positronium with zero total 
momentum I — and the large principal quantum number n. Its annihilation rate into 
photons to the leading order in QED couphng reads 

n + 1 

ToPsn^ST = 7 , (46) 

thus numerically, 

n + 1 

roPs„^37 - 7 X 10^-^ X ^ . (47) 

At large n it scales as oc n"^ and becomes much smaller than the decay rate of the or- 
thopositronium ground state. Meanwhile, the fastest transition from the level ns is to the 
level 2p, in which the rate is 

3 



^ ns-^2p 



4 








3 


[ 4 







x\d\\ (48) 



where m is the electron mass, and the squared matrix element of the dipole moment [42] is 
given by 

2_ 4 2i-V(n-2)2"-6 
' ' ~ m^a 3(n + 2)2»+6 " 

At a large n it approaches 

,,,2 4 2i5gxp(_8) 



ni?a 



3n3 



and thus the transition rate (48) asymptotes to 

The estimates in (47) and (49) explain what happens to the Rydberg positroniums. Low 
excited levels of n < 20 decay into 2p state, which quickly (r2p^is = 3.1 x 10^ s~^) decays 
further into Is, where finally positronium annihilates. At a larger n the direct annihilation 
of excited states dominates. Note that the matrix element of the positronium oscillation (3) 
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also depends on n. The estimate of the total oscillation probability of the excited positronium 
follows from the formula in (12) upon rescaling 

. .n+1 ^ n+1 

^^^^^ 7^r„p.„^3. = 7^. 

Hence, at large n the invisible decay branching ratio of the Rydberg positronium coincides 
with that of the ground state. 

Note that for the model with small vacuum splitting between positronium and its twin 
one can think of oscillations between states of different principal numbers n and n'. Then 
both E and E', and 7 in two diagonal entries depend on these numbers and the formula for 
the disappearance branching ratio can be generalized to this case. 



5 Conclusions and open problems 

In this paper we have presented the complete analysis of orthopositronium oscillations into 
its twin within the Mirror World models. We took into account the relevant effects due to 
the possible scatterings off walls and gas atoms, and the influence of external electric and 
magnetic fields. 

In a background-free case, the highest sensitivity to the mixing of a photon and parapho- 
ton responsible for the oscillations in positronium sector, is achieved in the pure vacuum with 
an infinitely large experimental volume. At a given value of mixing parameter the invisible 
branching ratio of positronium (disappearing via oscillation to its twin, which subsequently 
decays into paraphotons) decreases when applying electric or magnetic fields, adding gas to 
the cavity and decreasing its size (or equivalently increasing positronium velocities). We 
have calculated this decrease as a function of the relevant physical parameters, which can 
be used to get rid of possible systematics in the experiment, if evidence of a positronium 
disappearance is found. 

To perform numerical estimates, a computer code has been written, which allows us to 
account for a realistic geometry of the cavity, a realistic positronium spectrum (distribution 
over velocity), the details of scattering off the cavity walls and gas atoms, etc. To illustrate 
the code, the numerical results for the setup of positronium experiment [8] proposed within 
the AEgIS project have been presented. For the first time, the experiment will enter the 
phenomenologically interesting and cosmologically allowed region of the model parameter 
space. The code can be used for other setups and may be further improved by implementing 
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more details of positronium interactions with gas, walls, and magnetic and electric fields. The 
properly completed modification can be used in simulations of events in a real experiment. 

The presented numerical results for the experimental setup [8] can be improved with the 
account of high-velocity positroniums. Their presence in a small number seems to slightly 
diminish the sensitivity to the mixing parameter: faster scatterings off gas atoms and walls 
spoil oscillations to the Mirror World. A similar effect — decreasing sensitivity to invisible 
mode — is expected for a more accurate realistic treatment of positronium scatterings in- 
cluding induced positronium annihilation, small blind spots in the detector, etc. Thus, we 
propose our numerical estimates to be used to place an upper limit on the sensitivity of a 
given experiment in the mixing parameter. 

For the first time, we analyzed oscillations of excited positroniums — Rydberg positron- 
iums — into their similarly excited twins. In the vacuum the branching ratio of a high-level 
Rydberg positronium decay into nothing is the same as that of the ground state (at the 
same mixing). Rydberg positroniums will be available, e.(/., in the AEgIS experiment and 
can be used for better control over possible systematics. They might be of some interest in 
the models with a slightly violated Mirror symmetry, resulting in a small mass shift between 
electron and its twin. 
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